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Abstract- The influence of wave frequency variation on the anomalous cyclotron 
resonance u = u>Be + kvu interaction (ACRI) of energetic electrons with a ducted 
finite amplitude whistler-mode wave propagating through the so-called transient 
plasma layer (TPL) in the magnetosphere or in the ionosphere is studied both an- 
alytically and numerically. The anomalous cyclotron resonance interaction takes 
place in the case when the whistler-mode wave amplitude By/ is consistent with the 
gradient of magnetic field Bq. The region of phase space occupied by anomalously 
interacting energetic electrons (synchronous particles) is determined. The efficien- 
cies of both the pitch-angle scattering of resonant electrons and their transverse 
acceleration are studied and the efficiencies dependence on the magnitude and sign 
of the wave frequency drift is considered. It has been shown that in the case of 
ACRI occuring under conditions relevant to VLF-emission in the magnetosphere, 
the energy and pitch-angle changes of synchronous electrons may be enchanced by a 
factor 10 2 -=- 10 3 in comparison with ones for nonsynchronous resonant electrons. So 
the small in density group of synchronous particles may give significant contribution 
to a whistler- mode wave damping in TPL. 



1. INTRODUCTION 



The cyclotron resonance of charged particles with whistler-mode waves (CRI) is one 
among the basic mechanisms which govern a dynamics of these waves in the magneto- 
sphere and in the ionosphere. As an example, CRI is closely related to one of the most 
fascinating phenomena - so-called VLF-emissions triggering in the magnetosphere (see, 
for example, review papers of Molchanov, 1985; Omura et al, 1991; Helliwell, 
1993 and Rycroft, 1993 and papers of Dysthe, 1971; Karpman et al, 1974; Nunn, 
1974; Dowden et al, 1978; Matsumoto, 1979; Bell, 1984; Nunn, 1984). Moreover, 
CRI may give also the main contribution to VLF-waves damping and growth and caused 
by them the pitch-angle diffusion of magnetospheric plasma hot population to loss cone 
results to an anomalous precipitation of energetic electrons from the radiation belts into 
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the upper atmosphere. From practical point of view, registration of VLF-emissions in 
the magnetosphere and in the ionosphere and charged particle precipitation induced can 
be applied to various diagnostic problems, for instance, to global control of the geliogeo- 
physical enviroment, to the forecasting of radiowaves propagation in the Earth-ionosphere 
waveguide, to the satellite monitoring of natural crisis processes like typhoons and earth- 
quakes and so on. 

Specific feature of CRI in the magnetosphere is the key role of inhomogeneouty of the 
geomagnetic field B and the plasma density n e which determine the typical space scale of 
cyclotron resonance region, its location along the geomagnetic field line, the intensity and 
direction of energy transfer at the nonlinear stage of a wave-particle interaction (WPI). 
The estimates performed (e.g. Bell and Inan, 1981) have shown that under the typical 
magnetospheric conditions the energetic electron scattering by VLF-wave with significant 
changes of particle energy and pitch-angle occurs only if the wave amplitude becomes 
large enough because the typical space scale of CRI-region is usually much less than the 
inhomogeneity length of geomagnetic field Bq(S/Lb). Indeed, if oo and k are the whistler- 
mode wave frequency and wave vector respectively, Y = uiHe/w is the energetic electron 
dimensionless gyrofrequency, S - is the arc length along the geomagnetic field B (S/L B ) 
and Lb is the inhomogeneity length, then one can determine the small parameter of 
problem considered by the following formulae 

(l + «i) kL B 

where a\ = Vj_(Y — 1)/[(2Y + 1) v R ] - is the square of electron perpendicular velocity 
normalized on the typical value, v r = (u — uJuejjk is the resonance velocity. According 
to the linear theory of WPI, the typical space scale of cyclotron resonance region located 
outside the equatorial plane is of the order of Ir = L B 5 1 ^ 2 . At the equatorial plane, the 
geomagnetic field gradient becomes zero. So in the case of WPI located at the equato- 
rial plane, the space scale of cyclotron resonance region increases up to Ir = LbS 1 ^ 3 . 
Nevertheless, in both cases under conditions, typical for VLF-emissions generation in 
the magnetosphere, the CRI-region space scale Ir is about two order of magnitude less 
than the inhomogeneity length L B - Therefore, it is of considerable importance to study 
the possibility of sharp growth of CRI temporal duration, for instance, due to extent 
of the interaction region space scale over the substantial portion of the inhomogeneity 
length Lb- For the cyclotron resonance interaction of energetic electrons with the ducted 
whistler-mode wave of variable frequency in the equatorial plane vicinity, this problem 
was considered by Brinca (1981) and Bell and Inan (1981). It was shown that the 
wave frequency variation in the case of optimum frequency function allows to enchance 
significantly the interaction region space scale for the most stable trapped electrons. In 
the case of fixed-frequency wave this problem was studied by Erokhin (1995) and it 
was founded that the ACRI takes place in the stationary transient boundary layer under 
consistency of the whistler-mode wave amplitude B w with the magnetic field gradient. 
So the following condition must be fulfilled : B w /B ~ l/kYL B ■ This condition is 
in analogy with the long-lasting resonance condition considered by HELLIWELL (1967) 
and the second-order resonance one described by Nunn (1971) but in contrast to these 
articles, paper of Erokhin (1995) relates to so-called synchronous particles whose phase 
$ (its definition see below) is close to n/2 during their crossing of TPL. As a result in the 
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case of synchronous particles the cyclotron resonance interaction becomes a large-scale 
phenomena because the interaction region space scale is comparable with the magnetic 
field inhomogeneity length. In TPL it is observed the antidrift dynamics of synchronous 
particles and relative changes of their energy and pitch-angles are of the order of 100 
percents if the magnetic field variation is large enough 5B ~ B . Consequently, in the 
transient plasma layer the anomalous cyclotron resonance interaction of synchronous par- 
ticles with the whistler-mode wave takes place and one would expect ACRI to modify the 
wave damping. 

As both rising and falling tones are observed in the magnetosphere when trigger- 
ing VLF-emission, it is necessary to study the wave frequency variation influence on the 
anomalous cyclotron resonance interaction of energetic electrons with the ducted whistler- 
mode wave in the transient plasma layer. The present paper is devoted to solving this 
problem. Its solution allows to perform correct estimates of the wave damping for anoma- 
lous CRI in TPL. 

The paper structure is the following. The basic equations derivation and relations 
resulted are given in Section 2. The case of fixed wave frequency is described in Section 
3. Analytical and numerical results of studying the frequency sweeping influence on 
anomalous CRI in the stationary TPL are given in Section 4. Results obtained are 
discussed in Section 5. 



2. BASIC EQUATIONS. 



Let us consider the cyclotron resonance interaction between energetic electrons and 
the ducted wistler with a frequency uo < uoBe propagating along a weakly inhomogeneous 
both plasma and magnetic field B . As this interaction is localized in the vicinity of 
field line it is quite natural to use the curvilinear orthogonal coordinate system with the 
basic vector along the arc length S of the field line, normal and binormal to this line. It 
is convenient to introduce the dimentionless variables s = cuS/c, t' = ut, f3 = v/c, 
where v - is the velocity of electron. 

The equations of motion for the non-relativistic electrons mirroring in the magnetic 
field B took the standard form (see, for example, Dysthe,1971; Nunn, 1974): 

f = ^sin*-/^, *=fl. (1) 

Here flw = \e\Bw/Tn e cuj - is the dimensionless wistler amplitude, j3 p h = {uj/uj pe ){Y— 
1) 1/2 and P g = -2p R /Y - phase and group whistler mode wave velocities respectively, 
(3r = (1 — Y)f3 p h -resonance velocity, Y s = d s Y - magnetic field gradient, $ - the 
complement of the angle between the electrons perpendicular velocity v± and Bw- The 
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slow frequency sweeping rate of the whistler is taken into account by term 5uj/oo in the 
wave phase equation, and under wave amplitude diffusion neglecting it depends only on 



/ _ / f ds' 

Having in mind the VLF-emissions magnetospheric typical parameters one can assume 
the whistler-mode wave frequency change to be small enough (i.e. \5lj\ <C uj) during the 
particle crossing the resonance region. 

It is nesesary to pay attention to the following circumstance. As far as we are in- 
teresting in the WPI large-scale dynamics when the resonance region is global in a size 
comparable with the inhomogeneity length of static magnetic field B , in the equation 
(1) for phase $ we should keep the second term proportional to the whistler amplitude 
VL W . In the case of small-scale WPI analysis, it is usually neglected (see, for example, 
Matsumoto, 1979; Brinca, 1981). 

According to Erokhin (1995) and Erokhin et al. (1995), the magnetic field and the 
plasma density in the transient plasma layer are monotonous functions of the arc length 
S. So it is possible to put them in one to one correspondence ui pe = u pe (Y) , where ui pe 
is the electron langmuir frequency. To concretize the following calculations we are using 
the power function of the type u pe (Y) = uj pe (Yo)(Y/Yo) a , where Yq is the magnetic 
field at some point So inside TPL. For the power index a it is usually taken the value 
(j = (DE-model) or a = 0.5 (CL-model; see, for example, Bell and Inan, 1981). The 
spatial dependence of whistler amplitude flw can be given by the condition of energy 
flux conservation in the ray tube which crossection is inversely with B (see, for example, 
Molchanov, 1985). Therefore we obtain the following scaling of the cyclotron resonance 
and phase velocities as well as the whistler amplitude on the magnetic field Y : 

Pph = ^Jy) ' Pr{y) = ~^^T' (2) 

n w (Y) = Q*Y 1/2 Q(Y), Q(Y) 



where the following notations were used 

^=P ph (Y )Q 2 (Y ), 0,= 



" (y- 1)V4' 



Y 1/2 Q(Y ) 



For convenience of the subsequent analysis it is nesessary to transform the set of 
equations (1) to canonical form. Let us introduce the typical space scale of TPL as L t . 
In dimensionless variable s it is equal to = uL t /c = /3*/2f2* ^> 1 . Applying the scale 
transformations of variables f3\\ = f3*u, (3± = f3*v, s = s*£, t' = r/2fi* to the set of 
equations (1) and taking into account (2) we reduce (1) to the canonical form: 

2^ = Q(Y)Y 1 / 2 v sin^ - v 2 F(Y), ^ = u, 
dr dr 

dv Y 1 / 2 

2 d^ = Q(Y) [1 ~ UQiY)] $ + UVF{Y)l (3) 



4 




x [uQ\Y) + Y-l] + 



yl/2 



[1 -wQ 2 (r)]cos$ -x(l 



q(>> 



where x = (W > 1, « 9 = 2 (Y - 1)/Q 2 (F)F and F(Y) = Y 6 /Y is the magnetic 
field logarithmic gradient. Eqs. (3) form the basic set of equations for the subsequent 
analysis and the function determines the TPL spatial structure. It will be founded 
below from the condition of the synchronous particles existence. To describe briefly 
the TPL spatial structure in the case of constant wave frequency one puts 5uo = in 
Eqs. (3). According to the paper by EROKHIN et al. (1995), the anomalous cyclotron 
resonance interaction of energetic electrons with the ducted whistler-mode wave takes 
place for the group of synchronous particles defined by the following conditions: 1) the 
phase $ along the synchronous particle path is constant and is equal to $ s = n/2; 2) the 
parallel velocity of synchronous particle u s is equal to the cyclotron resonance velocity 
ur = — (Y — 1)/Q 2 {Y) . Putting in (3) u = ur, $ = 7r/2 and Sou = we obtain the first 
order nonlinear equation for the magnetic field profile in TPL: 



where /i — 1 + a/2, R(Y) = u R (Y) and A is a positive parameter determining the 
one-parameter set of the magnetic field profiles Y(A, £), for which the anomalous CRI of 
energetic particles with the small amplitude whistler-mode wave may occur at the entire 
transient plasma layer. 

Let the power index a be in the range < er < 3/2 . The solution of equation (4) 
exists if and only if the function F(A, £) is in the region (1, Y m (A)) , where Y m (A) is the 
single root of equation A = R(Y m ) , monotonously increasing under parameter A growth. 

To determine Y(£) we suppose the magnetic field in TPL to be varying in the range 
Y\ < Y(£) < Y2 , where Y\ > 1 and Y2 < Y m (A) . Then the magnetic field profile Y(£) 
can be obtained by inversion of the following monotonous function 



Fig. la depicts the part of profile (5) in the case when the magnetic field is varying in 
the range 2 < Y(£) < 4 , the power index corresponds to DE-model (a = 0.5) and the 
parameter A takes the following values: 6.755, 18 and 40. For the given values of param- 
eters a and A possible maximum magnitudes of the magnetic field Y m are respectively: 
4.001, 5.674 and 7.755. 

Fig. lb depicts the profile (5) with magnetic field variation in the range 2 < Y(£) < 4 
with A = 40 and different values of power index a. According to Fig. lb for the given 
value of parameter A, the TPL width is larger in the case of DE-model. 

For the transient plasma layer with magnetic field variation in the range Y 1 <Y(g) < 
Y 2 , the layer width is determined by the following expression: 



dY_ 



YF(Y) 



Y »(Y - 1)V 4 [A - RjY)} 1 / 2 
A + 2(1 - a + a/Y)R(Y) 



(4) 




Y [A + 2(1 - o + -)R{x)\ dx 
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Fig.la. Magnetic field profile V(0 i n the transient plasma layer for a = 0.5 and 
different values of parameter A: 1 - A = 6.755; 2 - A = 18; 3 - A = 40. TPL 
corresponds to the dispersed magnetic field jump 2 < Y(£) < 4 . 




Fig. lb. Magnetic field profile in the case 2 < < 4 for A = 40 and different 
values of power index a: 1 - a = 0.5 (CL-model); 2 - <j = (DE-model). 



or in the dimensional variable S we obtain L s = (cf3*/2u£l*) . Fig.lc illustrates 
the dependence of TPL-width 1% on parameter A, defined by formula (6) in the case of 
Y\ — 1, Y2 = ^m(A) for two values of the power index a. It can be seen again that 
the layer width is larger for DE-model of the plasma density. As far as for synchronous 
particles there is a well known integral of motion (Karpman et ai, 1974; Nunn, 1974) 



< + (1 - y>s = A » 

their energy S s , magnetic moment fi s and velocities components determined by the mag- 
netic field local strength 



u s (Y) 



YA - R(Y) 

' 2(y-i) ' 

(Y - I) 3 / 2 



W J_ 2F ~ 2(y-i) ' 



(7) 



.(y^T) 



.1/2 



(a - i?(y)) 
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Fig.lc. Dependence of TPL-thickness 1%(A) with magnetic field variation 1 < 
V(£) < V m (A) on parameter A in the cases: 1 - a = 0.5; 2 - a = 



The time Ar s , required for TPL-crossing by the synchronous particle, is equal to 

Y 2 



Ar.=/ 

Yi 



dY 



YF(Y) u g {Y) 



Let us consider again equations (3). To perform the subsequent calculations we assume 
the following model of wave frequency variation 8u/u = h(p) , where h(p) is the given 
function. The spatio-temporal dependence of function r) is governed by the following 
equation dp/dr — v{l — u/u g ) , where parameter v determines the frequency sweeping 
rate. The frequency sweeping rate can be characterized by a function T = (df/dt)/f 2 
which has the folowing scaling 



T ~ 4 • 10" 5 



df/dt 
1 kHz/sec 



5 kHz 

— T 



1 2 



orientated towards the magnetospheric VLF-emissions parameters. Functions T and h are 
related by dh/dp = yY jAmv . Therefore in the case of constant wave frequency sweeping 
rate, corresponding to the choice h{p) = p, parameter v is equal to v — x^/^n . So 
under the typical conditions of VLF-emissions in the magnetosphere, one has v < 1. 



3. ESSENTIAL FEATURES OF THE ANOMALOUS CYCLOTRON 
RESONANCE INTERACTION OF ENERGETIC ELECTRONS WITH THE 
FIXED FREQUENCY WHISTLER-MODE WAVE IN TPL. 



Before the studying the whistler frequency drift influence on the anomalous CRI in 
TPL, to clarify the following analysis it is nesessary to describe shortly the case of fixed 
frequency wave coresponding to the condition 5uj = in Eqs. (3). Let us assume that in 
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stationary TPL the magnetic field varies in the range Y 1 < Y(£) < Y 2 and the whistler- 
mode wave propagates in the direction of Y(£) growth. Therefore, the resonant particles 
are travelling in the opposite direction towards the wave. The maximum cyclotron reso- 
nance interaction in the transient plasma layer takes place for the group of synchronous 
electrons whose phases $ s are constant and equal 7r/2 and other parameters are defined 
by (7). The group of synchronous particles can also be defined in the different way. Let 
us introduce function J by the following expression: 

J = M 2 + (l-i), 2 . 

According to papers of Karpman et al. (1974) and Nunn (1974), for resonant par- 
ticles function J is the approximate integral of motion. Consequently, the group of syn- 
chronous particles can be defined by the folowing conditions: u s ~ ur, $ s ~ it/ 2 and 
J ~ A . Nonconserving part of the approximate integral of motion J can be easily esti- 
mated from (3). Performing the asymptotic integration, we obtain the modified integral 
of motion 

J = J+^cos$. 

Vy 

In turn, non conservation of the function I is caused by fastly oscillating terms of the 
order of 1/x 2 but the particle crossing of the local hyroresonance region is accompanied 
by the jump of function I proportional to l/x 3//2 . 

Let us turn back to formulae (7). According to (7) when crossing the TPL the syn- 
chronous particles increase their energy, magnetic moment, perpendicular velocity and 
pitch-angle a s = arctan (i; s /|u s |) = tan -1 (v s /\u s \) . At the same time, when moving in 
the direction of magnetic field decreasing, their parallel velocity goes down, i.e. in the 
transient plasma layer the antidrift dynamics of synchronous particles takes place. Earlier 
similar effects for electrons, trapped by the whistler-mode wave in the equatorial plane 
vicinity, were pointed out by Matsumoto (1979). 

Now let us assume that at TPL-entrance the velocities of incoming electrons coincides 
with the synchronous particle velocity, i.e. u(0) = u s (Y 2 ) , v(0) = v s (Y 2 ) , but there is a 
small deviation 6q in the initial phase $(0) = $ s + 0q . 

Numerical calculations and analytical estimates show us that there are some constants 
C(_) < < C( + ) so that for the initial phase detuning c/_\/x < 0q < ct + \/x during the 
resonance particle pass TPL, its phase is confined in the range < $ < n and its energy 
and pitch angle changes are close to the synchronous particle ones. If the initial phase 
detuning is large 9 > C(+)/% or 9 < C(_)/% , the duration of resonance interaction Atr 
becomes less than the time required for the resonance particle to cross the transient plasma 
layer. Under the growth of 6 Q the time Atr decreases as Atr ~ l/(x l^ol) 1 ^ 2 • So there 
is the following scaling of the resonance particle energy growth AS in dependence on the 
initial phase detuning 6 : AS ~ l/(x\%\) 1 ^ 2 ■ Fig. 2a depicts the chart of the resonance 
particle relative energy growth A = 10 2 A£/£ in dependence on the normalized initial 
phase detuning x@o for small 6q. The system (3) was integrated numerically for TPL with 
the magnetic field variation in the range 2 < Y(£) < 4 and the following parameters: 
a = 0.5, A = 10, v — . The initial data correspond to the synchronous particle i.e. 
Uo = u s (Y 2 ), v = v s (Y 2 ) , where Y 2 = 4 . According to Fig. 2a the maximum cyclotron 
resonance interaction of energetic electrons with the whistler mode of constant wave fre- 
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Fig. 2a. Dependence of synchronous particle energy gain A on normalized initial 
phase detuning x@o '■ 1 - X = 10 3 ; 2 - % = 10 4 . 

quency takes place for the small initial phase detuning \0 O \ ~ 1/x ■ It can be seen that 
charts are nonsymmetric on the phase detuning 8q and the resonant particle is gaining 
more energy in the case of positive 9$. In range x\@o\ > 10 2 the dependence A on the 
parameter x becomes appreciable. Similar effect is also observed for the pitch angle scat- 
tering of synchronous particles in TPL. Let Act represent the synchronous particle pitch 
angle change during the cyclotron resonznce interaction with wistler in TPL. We charac- 
terize the efficiency of electron pitch angle scattering during resonant WPI by function 
B = 10 2 (A«/o;o) , where ao is the electron pitch angle at the resonant region entrance. 
Notice that in theory of the VLF-emissions triggering in the magnetosphere the pitch an- 
gle scattering of energetic electrons by whistler-mode wave is usually characterized by the 
change of equatorial pitch angle a = sin -1 [(F e sin 2 a / F) 1 / 2 ] , where Y e is the magnetic 
field strength at the equatorial plane. The dependence of B on the normalized initial 
phase detuning x@o is shown in Fig. 2b for problem parameters as in Fig.2a. 




Fig. 2b. Dependence of pitch angle scattering efficiency B of synchronous particle 
on normalized initial phase detuning x@o '■ 1 ~ X = 10 3 ; 2 - x = 10 4 - 
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It can be seen that the pitch angle scattering efficiency B is larger approximately by 
factor 1.5 of the energetic efficiency A. From charts of Fig.2a and Fig.2b it follows that 
the cyclotron resonance interaction of synchronous particles, characterized by small 9 , 
with the wistler-mode wave in TPL is anomalously strong because of their energies and 
pitch angles relative changes are about 100%. Under 8q growth the efficiency of cyclotron 
resonance interaction falls. In the range of \9q\ ~ 1 the resonance region space scale is 
small enough. So with accuracy of a numerical coefficient of the order of unity one can 
use the linear estimate of A 



\A\ < 



sin2o Q 3/2 
Vx (A-y)V4 



A + 2 (1 - a + a/Y) R 
3 -2a + 2 a/Y 



1/2 



1 



(Y - I) 9 / 4 



(8) 



Taking x = 10 4 an d other parameters values correspondingly to Fig.2 we obtain from 
(8) that \A\ < 0.2 i.e. for large 6 the relative change of resonance particle energy is 
about 400 times less than the synchronous particle one. 

Results of numerical calculations of function A obtained by integrating system (3) in 
the case of large phase detunings # are given in Fig. 2c. The problem parameters are the 
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Fig.2c. Dependence of resonance particle energy gain A on initial phase 6q in the 
case of large 9q. The problem parameters are the following: x = 10 4 , v = 0, A = 
10, a = 0.5, Y = 4, u = u s (Y ) and v = v s (Y ) . 

same as for Fig.2a and x = 10 4 . Notice the narrow peak of the energetic efficiency A 
for small 6 corresponding to the synchronous particles. For comparison let us consider 
the cyclotron interaction of trapped electrons with the wistler-mode wave in TPL. The 
trapped particles perpendicular velocities at the TPL-entrance should be different the 
synchronous particle one v s (Y Q ). Introduce the normalized perpendicular velocity ae = 
v/v s (Y) and functions g(Y) and p(Y, ae) : 

MY) = - - fl < F » (9) 

33 

where condition < g < 1 takes place. From formula (9) it follows that in dependence 



on the variable ae the function p has minimum at ae = ae*(F) = J(l/g) — 1 which is 
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equal to minp = = 2 yjg(l — g) < 1 . The trapped particle dynamics corresponds to 
the motion of nonlinear oscillator in the potential well U(9) = Q%(pO — sin 9) , where 
&b = X^Vsiy) Y 1 / 2 Q 3 (Y) 1 4 is the bounce frequency square. The potential well exists 
only for p < 1 . In the case of0<g<l/2 the condition p < 1 is fulfilled for ae in the 
range 1 < ae < ae^ . If g(Y) is ranged as 0.5 < g < 1 , the trapped particle perpendicular 
velocity is less than the synchronous particle one, i.e. ae is in the interval ae^ < ae < 1 . 
The most long CRI of trapped particles with whistler in TPL takes place for the initial 
phase 9 = arccosp(F ) = ®b and when the trapped particle region has a 

maximum size in the phase plane (9, u) . 

The numerical calculations of the cyclotron resonance interaction of whistler-mode 
wave with the stably trapped particles with energy levels located closely to the potential 
well bottom were made, and the dependence of energetic efficiency A on the trapped 
particles perpendicular velocity ae = ae(Yo) was studied. Fig.3a depicts the chart of 
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Fig.3a. Dependence of energy gain of trapped particle, located closely to potential 
well bottom, on normalized perpendicular velocity aeo- Initial phase is equal to 
6» = cos" 1 p(ae ). 



A 




Fig. 3b. Dependence of resonance particle energy gain A on perpendicular velocity 
aeo and initial phase #o . The problem parameters are the following: \ = 10 4 > v = 
0, A = 10, a = 0.5, Y = 4, u = u s (Y ) and v = v s (Y ) . 
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function A dependence on aeo when v — 0, \ = 10 4 ) A = 10, o = 0.5 and Y = 4 . 
The trapped particle initial conditions were u = u s (Y ), 9 = #&(ae ) . The sharp 
peak of A(ae) near the point ae = 1 corresponds to the synchronous particles population. 
According to Fig. 3a, CRI- efficiency for the tapped particles is less by factor of few 
times than in the case of synchronous particles but it depends more smoothly on the 
perpendicular velocity aeo and the initial phase 6q. Therefore, the phase plane region 
occupied by trapped particles is substantially larger than one occupied by the synchronous 
particles. More evident it is demonstrated in Fig. 3b showing the efficiency A(a5 , 9 ) 
dependence on both initial parmeters. Please pay attention to high and narrow peak of 
A(aeo, #o) corresponding to the synchronous particles and the low but wide maximum in 
the case of trapped particles. 



4. ANALYSIS OF WAVE FREQUENCY DRIFT INFLUENCE ON RESONANT 
WPI IN TRANSIENT PLASMA LAYER. 

Consider the wave frequency drift influence on the anomalous cyclotron resonance 
interaction of energetic electron with a ducted whistler-mode wave propagating across 
the stationary transient plasma layer. Suppose the magnetic field Y(£) to be static and 
given by (5) and the wave frequency sweeping to be described by a function h{p). Below 
h{p) is taken as the linear function h{p) — p . So the frequency sweeping rate becomes 
constant. Contribution of the whistler frequency sweeping to the effective potential well 
U(p, 6) is characterized by a function r(r) 



d 



Y l l 2 Qv s dr 



p = v 



Y 



r + 



-(1 

Q 2 1 



dY 



Y 

2 ' 



(10) 



Y 



YF(Y) u g (Y) 



Notice that for the parametrs values typical for the VLF-emission generation in the 
magnetospher, the function r(r) magnitude is usually small. The wave frequency drift 
causes the following modification of function p determining the potential U 



p(ae, Y) = aeg + (1 



0)/ae 



Consider the case when condition 1 + r — g > is satisfied. Then p being a function 
of ae, has minimum at ae*(F) = [(1 + r — g) / g] 1 ^ 2 which is equal to minp = = 
2[g (1 + r - g)] 1/2 and p* < 1 if the condition r < (2g - l) 2 /4g is fulfilled. So in the 
case g — 1 < r < (2g — l) 2 / 4g there are trapped particles with perpendicular velocities 
ae in the range aei < ae < ae 2 , where ee lj2 are determined by 



aei 2 = 



1 ± ^{2g - I) 2 - Arg 

2g 



On the phase plane (9 T , 6) the phase trapping region has the maximum size at ae = ae* 
with its boundary defined by the following equation 



0.5 6 2 + VL 2 b (p*6» - sin 6) = Vl 2 b (p* arccos p*-Jl- p 
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For faller with the frequency sweeping rate satisfying to the condition r < g— 1 < , 
p(ae) becomes zero at the point se = ee c {Y) = [(g — l—r)/g} 1 ^ 2 . Therfore in this case there 
are always the trapped particles with the normalized perpendicular velocities ae in the 
range 883 < ae < 884 , where 

J(2g-l)*-4rg =F 1 
*» = 2 9 ' 

Consider the whistler frequency drift influence on the cyclotron resonance interaction 
of synchronous particle with initial data uq = u s (Yq), aeo = 1, Oq = ■ Introduce the 
specific phase value 9 V = (2 1^1) 1//2 defined by the wave frequency sweeping rate. 

In the case of riser r > with the frequency sweeping rate large enough the typical 
duration of synchronous particle CRI with the whistler may be estimated as Atr ~ 
4C3 / (xl r l ) 1/72 > where C3 is a constant of the order of unity. As a condition to suppress 
the anomalous CRI in TPL we take the following: Atr < 1/4 . Hence we obtain the 
restriction from below on the r-magnitude r > r c = 6.5 • 10 4 C|/x 2 and taking into 
account (5) it can be rewritten as a condition on a parameter v: v > v c , where v c is the 
critical value of v. Consequently, in the case of riser with the frequency sweeping rate large 
enough [y > u c ), substantial decreasing of the efficiency of synchronous particle cyclotron 
resonance interaction in TPL is occuring. In the range of initial phases \0q\ < 6 V , the 
energetic efficiency of CRI has the following scaling on parameter v: A ~ 1 / (x 2 ! 7 "!) 1 ^ 4 • 
If |6> | > 6 V , then this scaling is changed to one given in Section 3: A ~ 1 / (xl^ol) 1 ^ 2 • 

In the case of faller (r < 0), the synchronous particle becomes trapped one. Therefore, 
under the same but small initial phases 6q the duration of synchronous particle CRI 
for faller is substantially larger than one for riser. The range of synchronous particles 
perpendicular velocities is rather narrow <5ae = (ae — 1) <C 1 , so under the low frequency 
sweeping rate when |r| <C 1, the potential well U is shrinking U{9) m Vt 2 {6 [r + (2g — 
1) 5ae] + 9 3 /Q} . Introduce notations d= \r + (2g - 1) <5ae | and 6 d = (2d) 1/2 < 1 . If 
xOd 3> 1 the energetic efficiency can be scaled in the following way. When \0 O \ < 6d we 
obtain A ~ 1 / (x^d) 1/2 ■ In the range \8q\ > 64 the energetic efficiency is determined by 
the initial phase detuning A ~ 1 / (xl$o|) 1//2 • 

The computer simulations of the cyclotron resonance interaction efficiency under the 
whistler mode wave propagating across the TPL in the framework of equations (3) with 
h(p) = p were performed with taking into account the wave frequency sweeping. For 
a = 0.5, A = 10, x = 10 4 , Y = 4, 0„ = 0, u = u s (Y ) and v = v s (Y ) , the results 
are shown in Fig. 4 and Fig.5. Fig. 4a and Fig. 5a depict efficiencies A and B in the case 
of riser (0 < v < 10 3 ) . Fig. 4b and Fig. 5b are displaying them in the case of faller when 
< —v < 10~ 2 . For clarifying the efficiencies behaviour, the normalized variable xW] 1 ^ 2 
is used for the horisontal axis and plots of log A, \ogB versus log (xM 1 ^ 2 ) are presented. 
According to Fig. 4 and in correspondence with above developed theory function A has 
the plateau with the maximum value A m ps 75 in the range of small \u\ where the wave 
frequency drift influence is negligible. Outside this region the energetic efficiency A{y) 
falls approximately as A ~ 1 / (x 2 )^!) 1 ^ 4 under \v\- growth. For the positive v (riser) the 
break of plot takes place at the point ( x 2 W\) 1 ^ 2 ~ 2 . 
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Fig. 4a. Dependence of synchronous particle energy gain A on parameter x^ 1 ^ 2 i n 
the case of riser and x = 10 4 - 
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Fig. 4b. Dependence of synchronous particle energy gain A on parameter xlH 1 ^ 2 
the case of faller and x = 10 4 - 
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If v is negative the plateau is substantially wider (0 < (x 2 !^!) 1 ^ 2 < 7) . Besides at the 
point (x 2 !^!) 1 ^ 2 ~ 700 the local maximums of A and B arise with A ps 16 and B rs 51 
(see Fig. 5b). Comparing plots of functions A and 5 one can see that the pitch angle 
scattering efficiency B is substanially larger of the energetic efficiency A, in particular, 
in the plateau region one has B ps 100. Notice also that outside the plateau region the 
pitch angle scattering efficiency B isn't governed by the simple power-law decay on the 
variable xM 1 ^ 2 

The numerical calculations of the trapped particles dependence on the normalized 
perpendicular velocity ae and the wave frequency sweeping parameter v were made. Sim- 
ulations results are given in Fig.6 in the case of a — 0.5, A = 10, Y — 4 and x — 10 4 . 
The initial data of resonance particle correspond to the trapped electron located closely 
to the potential well bottom, i.e. Uq = u s (Y ), t> — &v s (Yq) and 6q = cos _1 (p(ae, Y )) . 
Consider the trapped particle energy gain in the case of its cyclotron resonance interaction 
with riser (see Fig. 6a). Notice the characteristic features of the normalized energy gain 
A. At first, there is the sharp peak for small v and <5ae = ae — 1 . Secondly, outside this 
peak with parameter v fixed, the energy gain A has maximum on the variable ae at some 
point ae m (z/) corresponding to condition p ps 1. In the range 1 < ae < ae m the trapped 
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Fig.5 Dependence of synchronous particle pitch angle scattering efficiency B on 
parameter xM 1 ^ 2 : a - r i ser { v > 0); b - faller [y < 0). 

particles are absent so the energy gain is very small. For ae > ae m the trapped particles 
exist and as ae growth the energy gain A is smoothly decaying. 

The energy gain of trapped particles in the case of their interaction with faller {y > 0) 
is plotted on Fig. 6b. It can be seen that in contrast to the riser case now the function 
A(ee, v) does not go down at the low Sse. Thus the strong cyclotron resonance interaction 
of trapped particles with faller is observed in the more wide range of parameter v = 10 2 v 
variation. For example, according to Fig. 6b in the case when aeo = 1 and v = —3, 
the energy gain of trapped particle with the initial phase detuning 6*0 ~ 0.46 is equal to 
A ps 65%. This particle crosses TPL in time At ~ 1.49 and the relative change of its pitch 
angle is B ~ 101% . Notice that for the wave frequency sweeping rate corresponding to 
T = 4 ■ 10 -5 and the dimensionless whistler amolitude = 1 / x — 10™ 4 the normalized 
parameter V is equal to v — —3.18 . 



5. CONCLUSION 



The principal conclusions of analysis performed above are the following: 
1. In the stationary transient plasma layer (5) the cyclotron resonance interaction of 
energetic electrons with the small amplitude ducted whistler mode wave of the fixed fre- 
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Fig.6 Dependence of relative energy change of trapped particle located close to 
potential well bottom on perpendicular velocity as and wave frequency sweeping 
parameter v = 10 2 f : a - riser [y > 0); b - faller {y < 0). 



quency is determined by three groups of resonance electrons: the synchronous particles, 
the trapped and nontrapped resonance ones. For the synchronous particles the cyclotron 
resonance interaction is global in nature because the resonanse region extends over the 
entire transient plasma layer. If the magnetic field B has the significant variation in TPL 
i.e. SBq ~ Bq, then the relative changes of synchronous particle energy and pitch angle 
in course of the cyclotron resonance interaction are anomalously large up to the order of 
100%. For the trapped particles the CRI-duration is less by a few times in comparison with 
the synchronous particles one but this population is more numerous than the synchronous 
particles group. Finally, in the case of nontrapped resonance particles under the typical 
magnetospheric conditions the cyclotron resonance interaction is weaker by the factor 
10 2 -T- 10 3 times in comparison with the synchronous particles CRI. In turn on the phase 
plane the nontrapped resonance electrons population is the most numerous. So it can 
be expected that in some range of the problem parameters all these groups of resonant 
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particles may give comparable contributions to the whistler damping (growth) rate in 
TPL. However detailed calculations of the whistler damping (growth) rate is a subject of 
a special consideration. 

2. In the cyclotron resonance interaction of energetic electrons with the whistler mode 
wave of variable frequency in the stationary TPL, for synchronous particles with the 
initial data uq = u s (Yq), vq ~ v s (Yq) and $ ~ 7r/2 there is the critical value of wave 
frequency sweeping rate such that above it [y > z/ ), the typical time scale of CRI decreases 
proportionally to 1 / (x 2 ^!) 1 / 4 where x 2 \ v \ > 1- There is a plateau in the dependence of 
CRI-efficiency on parameter x 2 \ v \ ( see Fig. 4 and Fig. 5). The plateau width is larger in 
the case of faller WPI. 

3. In the case of riser with moderate frequency sweeping rate corresponding to v < 1, 
the anomalous CRI in the transient plasma layer takes place for the trapped particles with 
perpendicular velocity in the range ae > 1. For the given v the maximum of energetic 
efficiency A is observed at ae suiting to the condition of the potential well being shrinked 
(p(ae, v) ~ 1) and function A decreases smoothly under ae-growth. 

4. In the case of energetic particles CRI with a faller in TPL, the trapped particles 
anomalous CRI region on the perpendicular velocity axis extends to the range a? < 1, 
where the maximum CRI-efficiency is observed and it corresponds to p(ae, v) ~ 1. Under 
parameters range considered in this paper the efficiency of trapped particles CRI with a 
faller is larger approximately two times than one in the case of riser. 

In addition it is nesessary to note the following. Above we studied the cyclotron reso- 
nance interaction of energetic electrons with the whistler mode wave of variable frequency, 
propagating across the stationary transient plasma layer. It was proved that there are 
the synchronous particles groups whose CRI in TPL is of the global nature if the whistler 
frequency sweeping rate is low enough. In the opposite case the synchronous particles 
existence may be related to the following factors: a) the magnetic field B nonstation- 
arity; b) the specific modulation of whistler frequency in analogy with one considered in 
paper of BRINCA (1981) and in paper of BELL and INAN (1981). In both cases a) and 
b) there is a question about the temporal duration of ACRI under a continuous injection 
of energetic particles into TPL. Besides, in the case a) it is also necessary to consider the 
source of magnetic field Bq(s, t) nonstationarity and to take into account the influence of 
electric field E driven by the time- variable magnetic field B . 
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